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ABSTRACT 

Recent theoretical and observational studies on the assembly of early-type galaxies (ETGs) 
point towards an inside-out growth of their stellar mass characterized by extended low mass- 
density halos grown around compact and dense cores. Models can form ETGs at high-z as 
compact spheroids that then grow in size through dry minor mergers. Dry mergers would 
affect mainly the outskirts of the galaxy enlarging the size (i.e. the effective radius) keeping the 
inner parts and the total stellar mass nearly unchanged. Hence, the central stellar mass density 
i-^ ■ will not change with time in contrast to the stellar mass density within the effective radius 

which should decrease with time as the effective radius increases. Some previous observations 
are interpreted as supporting inside-out growth as the central stellar mass density of high-z 
ETGs is found similar to that of local ETGs. In this paper we derive the central stellar mass 
density within a fixed radius and the effective stellar mass density within the effective radius 
for a complete sample of 34 ETGs morphologically selected at 0.9 < Zspec < 2 and compare 
them with those derived for a sample of ~ 900 local ETGs in the same mass range. We find 
that the central stellar mass density of high-z ETGs spans just an order of magnitude and is 
similar to that of local ETGs, as found in previous studies. However, we find that the effective 
stellar mass density of high-z ETGs spans three orders of magnitude, exactly as the local 
ETGs and that it is similar to the effective stellar mass density of local ETGs showing that 
it has not changed since z ~ 1.5, in the last 9-10 Gyr Thus, the wide spread of the effective 
stellar mass density observed up to z ~ 1.5 must originate earlier, at z > 2. Furthermore, 
we show that the small scatter of the central mass density of ETGs compared with the large 
scatter of the effective mass density is simply a peculiar feature of the Sersic profile and hence 
is independent of redshift and of any assembly history experienced by galaxies. Thus, it has no 
connection with the possible inside-out growth of ETGs. Finally, we show a tight correlation 
between the central stellar mass density and the total stellar mass of ETGs in the sense that the 
central mass density increases with mass as Ai~^'^. This implies that the fraction of the central 
stellar mass of ETGs decreases with the mass of the galaxy. These correlations are valid for 
the whole population of ETGs considered independently of their redshift suggesting that they 
originate in the early-phases of their formation. 

Key words: galaxies: evolution; galaxies: elUptical and lenticular, cD; galaxies: formation; 
galaxies: high redshift 



1 INTRODUCTION 

The formation and the mass assembly history of eaiiy-type galax- 
ies (ETGs; ellipticals and lenticulars) are central issues in the gen- 
eral topic of galaxy formation, issues that are not fully understood. 
Much effort has been devoted to understanding the formation and 
the evolution of ETGs in the context of hierarchical models of 
structure formation (e.g. Khochfar and Burkert 2003; De Lucia et 
al. 2006; Hopkins et al. 2010). In these models the formation of 
ETGs is directly linked to mergers events. The old stellar popula- 
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tions observed in almost all local ETGs suggest that most of the 
stellar mass formed a very long time ago and not during recent 
episodes of star formation (e.g. Thomas et al. 2005; Renzini 2006). 
This led to conclude that mergers that formed ETGs either occurred 
at very high redshift, z > 4-5, concurrently with or even triggering 
large episodes of star formation (e.g. Naab et al. 2007; Khochfar S. 
and Silk I. 2006a) or were "dry", that is the coalescence of pre- 
existing old stellar systems (e.g. De Lucia et al. 2006). In fact, the 
current and widely accepted view involves a combination of these 
two merger types acting at different time. Indeed, ETGs are thought 
to form at high redshifts as compact spheroids resulting from a 
main gas-rich merger event in which most of the stellar mass of the 
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galaxy forms. They should then grow in size mostly by the addition 
a low stellar mass density envelope through subsequent dry minor 
mergers at later times (e.g. Hopkins et al. 2009; Naab et al. 2009; 
Bezanson et al. 2009). Naab et al. (2007) found that their simulated 
galaxies start forming stars at z ^ 3 - 5 concurrently with the in- 
tense phase of dissipative merging. In about 1.5 Gyr (at z ~ 2.5) 
almost 80 per cent of their final stellar mass is formed in the dis- 
sipative central starburst (see also Sommer-Larsen and Toft 2010) 
justifying the old stellar populations observed in local ETGs. Then, 
the resulting compact and massive spheroid grows in size through 
dry minor mergers whose main elfect would be that of re-arranging 
stars mostly those far from the galaxy center, efficiently enlarging 
the size but keeping the inner parts and the total stellar mass of the 
galaxy nearly unchanged (e.g. Naab et al. 2009). 

This scenario was triggered by the fact that the first observed 
high-z ETGs were characterized by an effective radius smaller than 
the mean radius of local counterparts with comparable mass, re- 
quiring a growth of their size across the time to match the size of 
the local ETGs. The above scenario, often called inside-out growth, 
qualitatively reproduces this possible size evolution of ETGs dis- 
cussed and debated in many studies (Daddi et al. 2005; Trujillo et 
al. 2006; Longhetti et al. 2007; Trujillo et al. 2007; Cimatti et al. 
2008; van der wel et al. 2008; van Dokkum et al. 2008; McGrath et 
al. 2008; Damjanov et al. 2009; Saracco et al. 2009; van der Wel et 
al. 2009; Trujillo et al. 201 1; Damjanov et al. 201 1; van de Sande et 
al. 2011). However, it requires a considerable fine tuning to repro- 
duce both the local scaling relations and their scatter which is sig- 
nificantly smaller than the one produced by minor mergers (Nipoti 
et al. 2009; Nair et al. 201 1). In addition to this argument, there are 
some observations that provide evidence that contradicts this stel- 
lar mass growth and the more general hierarchical scheme. Indeed, 
recent studies focused on brightest cluster galaxies at 1.2 < z < 1.5 
show that they were almost fully assembled by this time having 
already grown to more than 90 per cent of their final stellar mass 
at that redshift (Colfins et al. 2009; Stott et al. 2010). In agree- 
ment with these results, studies on the evolution of the galaxy stel- 
lar mass function show that massive ETGs were already assembled 
at z > 1 and that they have not grown further at lower redshift 
(e.g. Pozzetti et al. 2010). These results suggest that merging has 
had little effect on the stellar mass growth of ETGs in the last 9 
Gyr contrary to hierarchical models that predict a protracted mass 
build-up over Hubble time. Furthermore, in the last couple of years 
evidence of a large number of ETGs at z ~ 1.5 with sizes similar 
to those of local ones with comparable mass, here termed normal 
ETGs, has been accumulated showing that most of the high-z early- 
types were not more compact than local ones (Saracco et al. 2009; 
2010; Mancini et al. 2010). 

One of the differences between those works finding almost 
exclusively compact ETGs at high-z and those finding a majority 
of normal ETGs is the definition of early-type galaxy and, con- 
sequently, the selection criterion used to construct the samples of 
ETGs, a criterion that also determines the completeness of the sam- 
ples, possible selection bias and the redshift range covered. The 
studies showing that the population of ETGs at high-z is composed 
of compact galaxies are based on samples of galaxies selected ac- 
cording to their colors and/or on samples of passive galaxies se- 
lected according to their low value of the specific star formation 
rate, a quantity derived from the best fits of their spectral energy 
distribution with stellar population synthesis models. These selec- 
tion criteria allow the selectio of ETGs, according to the above def- 
initions of ETG, even beyond z ~ 2.5 (e.g. Buitrago et al. 2008; 
Cameron et al. 2011; Cassata et al. 2011). The studies showing a 



predominance of normal ETGs at high-z are based on morphologi- 
cally selected samples and, for this reason, are usually limited to a 
slightly lower redshift (z ~ 2) than the former (e.g. Saracco et al. 
2010; Mancini et al. 2010). We will discuss the implications that 
these different selection criteria of ETGs may have on the results in 
later sections of this paper. 

Concurrently with the evidence of a large number of normal 
ETGs at z ~ 1.5, evidence of the presence of a large fraction of 
compact ETGs in the local Universe similar to the high-z ones has 
emerged (Valentinuzzi et al. 2010a; 2010b). In fact, the number 
density of compact ETGs morphologically selected at z ~ 1.5 is 
found to be consistent with the number density of ETGs selected 
in the same mass range in the local Universe (Saracco et al. 2010) 
showing that size evolution, if it exist, cannot affect the majority 
of the high-z ETGs, in agreement with what is found from other 
analysis (e.g. Newman et al. 2010; Stott et al. 2011). These recent 
results indirectly cast some doubt on the inside-out growth as a 
general mechanism of accretion of stellar mass and of growth in 
size. On the other hand, other studies focusing on the central and 
the effective stellar mass density of high-z ETGs have shown that 
the former is almost constant and similar to that observed in the 
local ETGs (e.g. Bezanson et al. 2009; van Dokkum et al. 2010; 
Tiret et al. 201 1) providing one of the strongest arguments in favor 
of the hypothesized inside-out growth of ETGs throughout their 
life. 

In this paper we focus on the central and the effective stel- 
lar mass density of ETGs and on their relations with other fun- 
damental quantities. We use both a complete sample of ETGs at 
0.9 < Zspec < 2 and a local sample of ETGs selected according to 
the same criteria as used for the high-z sample in order to check 
for possible differences. The aim of this analysis is twofold. First, 
we want to asses whether the central and the effective stellar mass 
densities of ETGs are linked to the evolutionary path followed by 
ETGs and hence whether they can be used to probe their assembly 
history or whether they are independent on it. Secondly, we want 
to asses whether inside-out growth is a viable way to justify the 
properties observed in ETGs both at high-z and in the local uni- 
verse examining both the evolution they underwent at z < 1.5 and 
the early phases they possibly experienced at z > 2. In Sec. 2 we 
describe the data. In Sec. 3 we define the central and the effec- 
tive stellar mass densities probing the inside-out growth scenario 
at z < 1.5. In Sec. 4 we use the results obtained to constrain the 
general spheroids formation scheme at very early epochs. In Sec. 5 
we define and study a scaling relation involving the central and the 
total stellar mass, scaling relation followed by ETGs independently 
of their redshift. Finally, in Sec. 6, we summarize the results and 
present our conclusions. Throughout this paper we use a standard 
cosmology with Hq = 70 Km s"' Mpc"', £1„ = 0.3 and Qa = 0.7. 
All the magnitudes are in the Vega system, unless otherwise speci- 
fied. 



2 THE DATA SET 

The sample of ETGs we used in our analysis is composed of 
34 galaxies morphologically selected from the southern field of 
the Great Observatories Origins Deep Survey (GOODS-South; Gi- 
avalisco et al. 2004) described in Saracco et al. (2010). The sam- 
ple was constructed by first selecting all galaxies brighter than 
K^20.2 over the ~ 143 arcmin- of the GOODS-South field and 
then by removing all the galaxies with measured spectroscopic red- 
shift Zspec < 0.9 and those with irregular or disk-like morphology. 
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This first step of the morphological classification was undertaken 
through a visual inspection of the galaxies carried out indepen- 
dently by two of us on the ACS images in the F850LP band. Only 
galaxies for which the two independent visual classifications agree 
were removed. Then, on the basis of the best-fitting procedure to 
the observed profile described below, galaxies having a Sersic in- 
dex n < 2 01 clear irregular residuals resulting from the fit were 
removed. Out of the 38 early-type galaxies thus selected 34 had 
measured spectroscopic redshift 0.9 < Zspec < 1-92 leading to ~ 90 
per cent spectroscopic completeness. 

We note that at least eight galaxies of our sample show weak 
on emission line (see Tab. 1). The lack of X-ray emission asso- 
ciated to these ETGs suggests that the emission is most proba- 
bly associated with star formation rather than with AGN. Hence, 
these ETGs are not passive. Given the S/N (< 10) of the spectra 
at Arest -3500A-4000A and the continuum emission at these wave- 
lengths, we estimate that only emission stronger than ~ 10""* erg 
cm"^ s"' could be detected. Thus, the star formation rate for these 
eight galaxies would be in the range 0.4-2 solar mass per year, 
according to their redshift, where we used the relation SFR(M0 
year-')=(1.4 ± 0.4) x lO"'*' L[OII](erg s"') (Kennicutt 1998). It is 
worth noting that for masses of a few times 10'" Mq this would 
result in a specific star formation rate -11 <log(SSFR)< -10 yr"', 
higher than the SSFR often used to define passive galaxies (e.g. 
log(SSFR)< -11 in Cassata et al. (2011)). We cannot exclude that 
the number of non-passive ETGs in our sample is higher than eight. 
Indeed, we note that the spectroscopic observations collected for 
these eight galaxies (di Serego Alighieri et al. 2005; Cimatti et al. 
2008; van der Wei et al. 2005) are deeper than those of the other 
GOODS spectroscopic observations. Consequently, for many of 
them,in particular for those at z > 1 .4, spectroscopic observations 
do not allow detection of Oil emission lines fluxes as faint as a few 
10"'* erg cm"- s"', corresponding to SFR>> 1 Mg yr"' at that red- 
shift. Hence, the presence of weak star formation would be hardly 
detectable. We note also that ETGs with colors bluer than those 
often used to select passive galaxies at z > 1 (e.g. (R-K)> 5, or 
(I-H)Ag > 2.1, Cameron et al. 2011) are present in our sample. All 
these ETGs, those non-passive and those bluer than the above color 
cuts, would not be included in a sample of ETGs selected according 
to their SSFR and/or to their colors. The resulting 34 ETGs of our 
sample are at redshift 0.9 < Zspec < 1-92 and have stellar masses 
in the range 0.6 x 10'" < M, < 3 x 10" Mq as derived by fitting 
their spectral energy distribution (SED) sampled by 14 photometric 
points in the wavelength range 0.3-8.0 //m (Saracco et al. 2010). 



3 THE CENTRAL AND EFFECTIVE STELLAR MASS 
DENSITIES OF ETGS 

3.1 Surface brightness profile fitting 

The effective radius R,, [kpc] (r^ [arcsec]) of our galaxies has been 
derived by fitting a Sersic profile 
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to the observed light profile in extremely deep (40-100 ks) HST 
ACS-F850LP images. The point spread function (PSF) to be 
convolved with the Sersic profile has been constructed for each 
galaxy by averaging the profile of some unsaturated stars. The two- 
dimensional fitting has been performed using Galfit software (v. 
2.0.3, Peng et al. 2002). The fitting to the galaxy profiles provided 
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Figure 1. The effective radius of the 15 galaxies of our sample in the redshift 
range 0.96 < z < l.l.derived in the ACS-F850LP band (/i„.„ ~ 0.42 fim) is 
compared to the effective radius derived in the ACS-F606W band (A„st ~ 
0.3 //m). 

us with the axial ratio b/a and the semi-major axis of the pro- 
jected elliptical isophote containing half of the total light which is 
related to the circularized (averaged) effective radius by the rela- 
tion = \fbja. The effective radii R^, we obtained through the 
two-dimensional fitting procedure are in the range 0.4-9 kpc. 

A concern over the effective radii thus derived could arise 
in relation to the different rest-frame wavelengths sampled by the 
F850LP filter given the wide redshift range of our sample. In partic- 
ular, while at z ~ 1 the F850LP filter samples the rest-frame close to 
the B band (~ 0.42 fim), at z > 1.5 this filter samples the rest-frame 
at ~ 0.3 ytim, close to the U and NUV bands. However, the emission 
in the B and U bands of a galaxy is dominated by the light coming 
from the same young stellar component, in contrast to the emission 
at A„st > 0.65-0.7 jim (R band) which is dominated by the old stel- 
lar component. Hence we do not expect a dependence of the size 
of our galaxies on the different rest-frame wavelength sampled, i.e. 
on the B-band or U-band according to the redshift of the galaxy. In 
order to exclude the presence of this dependence we selected the 
15 galaxies of our sample in the redshift range 0.96 < z < 1.1 and 
we compared their effective radius measured in the F850LP band, 
sampling /i„„ ^ 0.42 fim, with the effective radius measured in the 
F606W band, sampling /l„„ ^0.3 yum. The comparison is shown in 
Fig. 1. It can be seen that no dependence on wavelength is present 
as expected. 

3.2 Stellar mass density estimates 

In Fig. 2, the effective radius R^, of the 34 ETGs of our sample 
(filled symbols) is plotted as a function of their stellar mass M,. 
The solid line is the local size-mass (SM) relation of Shen et al. 
(2003), while the dotted lines mark one sigma scatter. According 
to the definition already used in Saracco et al. (2010), cyan points 
represent the ETGs falling within one sigma from the local SM 
relation, here termed normal, while magenta points are compact 
ETGs, those falling at least 1 sigma below the relation. Crosses 
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Figure 2. Upper panel: the effective radius of the 34 ETGs at 0.9 < 
ZsiMc < 2 (filled symbols) and of the ~ 900 ETGs selected from the WINGS 
survey (crosses) is plotted as a function of their stellar mass M, . The solid 
line is the local size-mass (SM) relation by Shen et al. (2003) with its 1 
sigma scatter (dotted lines). Cyan points are those high-z ETGs (normal) 
falling within one sigma from the local SM. Magenta points mark high- 
z ETGs more compact falling at least 1 sigma below the relation. Lower 
panel: the compactness, defined as the ratio between the effective radius Rf 
of the galaxy and the effective radius R(..7=o of an equal mass galaxy at z=0 
as derived by the local SM relation, is plotted as a function of the stellar 
mass. Symbols are as in the upper panel. 



represent ~ 900 ETGs we selected at z < 0.05 from the Wide- 
field Nearby Galaxy Cluster Survey (WINGS; Fasano et al. 2006; 
Valentinuzzi et al.2010a). We point out that the local SM by Shen 
et al. is taken as reference to define normal and compact galaxies. 
It acts in the same way on any sample of ETGs considered. For 
instance, we used this relation to define compact ETGs both in our 
high-z sample and in the WINGS sample to compare in a consistent 
way their co-moving number density (Saracco et al. 2010). The 
local sample of ETGs has been selected from the WINGS sample in 
the same stellar mass range and on the basis of their morphology, as 
done for our high-z sample. In particular we selected ellipticals (E) 
and lenticulars (SO) according to the morphological classification 
of the WINGS survey (Fasano et al. 2012). It is worth noting that 
WINGS effective radii have been measured on the V-band images 
(^rest - 0.5 yum; Pignatelli et al. 2006), very close to our /i„„ ~ 0.42 
pm). In the lower panel of Fig. 2, the compactness, defined as the 
ratio between the effective radius of the galaxy and the effective 
radius Re,z=o of an equal mass galaxy at z=0 as derived by the local 
SM relation, is plotted as a function of the stellar mass for both the 
samples. 

Integrating the best-fitting Sersic profile of eq. (1) over a pro- 
jected area A = nR\ we derived for each galaxy the luminosity L\ 
interior the radius R, = 1 kpc 



L\ = 2nIiR'^n-g^y{2n, x) 



(2) 



where n is the Sersic's index, x = b„(R] /Re)'^" and y(2n,x) is 
the incomplete gamma function. The choice of 1 kpc radius (2 
kpc diameter aperture) within which to calculate the central stellar 



mass density has been made in consideration of the resolution of 
the ACS-F850LP images characterized by a FWHM=: 0.15 arcsec 
corresponding to about 1.2 kpc at z ~ 1. However, it is important 
to note that the stellar mass densities are derived from the intrinsic 
galaxy profile, the Sersic profile. Thus, PSF does not affect in any 
way the estimate of the stellar mass densities. PSF plays a role only 
in the 2-dimensinal fitting to the observed profile as it is convolved 
with the Sersic profile. By replacing in eq. (2) 7(2n, x) with the 
complete gamma function r(2n) we obtained the total luminosity 
L,o, (Ciotti 1991) and the fraction of luminosity interior to Ri is 
then 



Li_ _ y(2n, x) 
L,o, T(2n) 



(3) 



In our computation we assumed the analytic expression b„ = 
I.9992« - 0.3271 (Capaccioli 1989) to approximate the value of 
b„. Assuming that the light profile traces the stellar mass profile 
and that the mass-to-light ratio (M/L) is radially constant across 
the galaxy, the stellar mass Mi interior to 1 kpc is given by Mi = 
L|/L,o, X M,, where M, is the total stellar mass of the galaxy ob- 
tained from the SED fitting. Because the radius Ri is fixed, the stel- 
lar mass density interior Ri is pi = cM\, where c = (4/3;r/?j)"' = 
0.239 [kpc"3]. Hence, pi and M\ follows the same behavior of the 
other quantities. Under the same assumptions we also derived the 
effective stellar mass density p^ = 0.5A1,/(4/3;r/?^), that is the den- 
sity interior the efi"ective radius. The whole set of parameters used 
and derived in this work for the sample of 34 ETGs is summarized 
in Tab. 1 . 

In fact, the mass-to-light ratio could not be radially constant 
for all the galaxies for the presence of color gradients in some 
high-z ETGs (e.g. Gargiulo et al. 2011, 2012; Guo et al. 2011). 
We investigated whether this can affect our estimate of stellar mass 
densities and our results. (U-R)„„ color gradients in high-z ETGs, 
when present, are negative indicating a redder color toward the 
center of the galaxy. Gargiulo et al. (2012), using synthetic stel- 
lar population model fitting to different galaxy regions, find that 
the main driver of such color gradients is the age of the stellar pop- 
ulations, older toward the center than in the outskirts with a possi- 
bly higher metallicity. Older ages imply larger M/L values, higher 
metallicity lower M/L values. Thus the two parameters act in an 
opposite way regarding M/L canceling out each other, at least par- 
tially. Anyway, let us assume that only the age varies radially thus 
producing the maximum M/L radial variation. According to this, 
M/L should be larger toward the center with respect to the M/L 
we derived from the SED fitting averaged over the whole galaxy. 
Consequently, both the central stellar mass density and the effec- 
tive mass density should be slightly higher. On the basis of the age 
gradients estimated by Gargiulo et al. (2012) we derived that the 
possible M/L variation is in the range 1-1.3. Thus, the central and 
the effective mass densities could be larger at most by this factor. 
In particular, for those galaxies whose effective radius is compara- 
ble to Ikpc (~ 12 galaxies in our sample) the effect will be exactly 
the same for both the stellar mass densities, central and effective. 
For the remaining galaxies whose effective radius R^ >> 1 kpc the 
effect will be slightly larger for the central stellar mass density than 
for the effective density. Gargiulo et al. (2012) detect color gradi- 
ents in 7 out of the 14 ETGs studied and find that color gradients 
are not related to the compactness of the galaxies. Thus, we expect 
that about 50 per cent of our sample show color gradients, inde- 
pendently of their compactness. It follows that, in the hypothesis of 
maximum M/L radial variation (i.e. of pure age variation) the dis- 
tributions of the central and effective stellar mass densities could 
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Figure 3. The distributions of tlie stellar mass density within the effective radius (left panel) and of the central stellar mass density pi within 1 kpc radius 
(right panel) are shown for the whole sample of 34 ETGs at 0.9 < z < 2 (soHd histograms). The shaded histograms represent the distributions of normal ETGs 
falling within one sigma from the local size-mass (SM) relation (cyan) and of ETGs more compact falling at least 1 sigma below the relation (magenta). 



be oifset by a factor in the range 1-1.3. Thus, the assumption that 
M/L is radially constant for all the galaxies does not introduce any 
significant systematic in our analysis even in the case that it is not 
strictly true. 

In Fig. 3 the distributions of the stellar mass density within 
the effective radius (left panel) and of pi within Ikpc radius (right 
panel) are shown. The black solid-line histogram represents the dis- 
tribution of the whole sample of 34 ETGs while the shaded his- 
tograms show the distribution of normal ETGs (cyan) and of com- 
pact ETGs (magenta). The distribution of the effective stellar mass 
density of Fig. 3 (left panel) varies by about three orders of mag- 
nitude, from ~ 10* Mo kpc"^ to ~ 10" Mq kpc"' with a median 
value = 7 X 10** Mq kpc"', reflecting the large variation (a fac- 
tor ~8) of for ETGs of the same stellar mass visible in Fig. 2. In 
contrast, the distribution of the central mass density pi shown in the 
right panel of Fig. 3 spans just an order of magnitude centered at 
Pi ^5x10' Mq kpc"'. The values of pi are actually within a factor 
2 from the median value as noted by other authors (e.g. Bezanson 
et al. 2009; Tiret et al. 2011). Analogously, while the median value 
p, of normal ETGs (3.5 x 10* Mq kpc"') differs by a factor 10 from 
the median value p^ of compact ETGs (3.9 x 10' Mq kpc"') in the 
case of pi the median values differ by less than a factor 2 (4 x 10' 
vs 7 X 10' Mq kpc"'). This information is summarized in Fig. 4 
where the central mass density pi of the 34 ETGs of our sample 
is plotted versus the effective mass density p, (filled circles). Our 
data are best-fitted by the relation pi oc p^ in agreement with 
Tiret et al. (2011) (to an increase by two orders of magnitude of 
P(, corresponds an increase by a factor 1.5 of pi) even if the cor- 
relation between the two quantities is not statistically significant 
being the probability of the Spearmen rank test P = 0.27. The 
small variation of the central mass density of ETGs with respect 
to the large variation of the effective stellar mass density is often 
used to advocate the inside-out growth of stellar mass experienced 



by ETGs through time (e. g. Bezanson et al. 2009; van Dokkum et 
al. 2010). In this scenario ETGs, after having assembled at high-z 
as compact/dense spheroids, gain a low stellar mass density enve- 
lope through subsequent dry minor mergers at later times growing 
efficiently in size and less in mass and reaching at z = the typi- 
cal size of local ETGs (e.g. Hopkins et al. 2009; Naab et al. 2009). 
In a plot like the one shown in Fig. 4, high-z ETGs would pop- 
ulate the right-hand side (high effective stellar mass densities) of 
the \Pe,Pt] plane and should migrate towards lower effective stellar 
mass density, the left-hand side, at lower redshift keeping almost 
constant the central density pi . In fact, some works based on in- 
complete sample of high-z ETGs pre-selected on the basis of their 
red color show this kind of segregation (e.g. Bezanson et al. 2009). 
On the other hand, it has been shown that this selection criterion is 
strongly biased toward compact ETGs (see Fig. 2 in Saracco et al. 
2010). To directly verify how color selection affects this plot, we 
selected the 10 reddest galaxies of our sample corresponding to a 
color (F850LP-K)^B > 2.4. These galaxies are marked by a black 
filled square in Fig. 4. It is evident that they preferentially populate 
the right-hand side of the plot and thus that a color selection pro- 
duce the segregation above. In contrast, it can be seen that the 34 
ETGs morphologically selected at z ~ 1.5 of our sample are not 
segregated in the right-hand side of the plot but rather they span 
three order of magnitude in effective stellar mass density in spite of 
their high redshift. Furthermore, in Fig. 5 the effective radius, the 
stellar mass, the central and effective stellar mass density of the 34 
galaxies of our sample are plotted as a function of their redshift. It 
can be seen that there are no significant correlation between these 
quantities and the redshift, as also confirmed by the Spearmen rank 
test. Hence, our analysis and our conclusions cannot be affected by 
any possible dependence on redshift of these quantities. Actually, 
Fig. 5 suggests also that there is no significant evolution of and 
Pp over the redshift range 0.9 < z < 2. Thus, because at z ~ 1.5 
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Figure 4. The stellar mass density within 1 kpc radius pi is plotted versus the effective stellar mass density p^. Large filled dots represent the ETGs of our 
sample at 0.9 < z < 2 (cyan noiTnal ETGs, magenta compact ETGs). Black filled squares superimposed to filled dots mark the 10 reddest galaxies of our 
sample, i.e. galaxies with color (F850LP-K)ab > 2.4. Black crosses represent the local ETGs selected from the WINGS sample. The open squai'es represent 
the 23 super-massive (0.3 - 5 X 10'~ Mq) ETGs selected by Bernardi et al. (2006; 2008) and studied by Tiret et al. (201 1). The red lines represent the locus 
occupied on the [pe.pi] plane by galaxies described by a Sersic profile with index n = 4, effective radii in the range 0.5 < < 10 kpc and stellar masses lO'" 
Mq (lower line), lO" Mq (middle line) and 4 X lO" M© (upper line) respectively. The dashed line is obtained for an index n = 2 and a mass 10' ' Mq. 



co-exist compact ETGs and normal ETGs, if we accept the frame- 
work of the inside-out accretion of stellar mass we must conclude 
that the latter must have increased their radius at earlier epochs, at 
z > 2, while the former will increase their size later, at lower red- 
shift. In the next section we compare this observational evidence 
with the inside-out growth accretion proposed to assemble ETGs. 

4 PROBING THE INSIDE-OUT GROWTH OF ETGS 

4.1 Are the central and the effective stellar mass densities 
related to the assembly history of ETGs ? 

Whether or not a fraction of the local ETGs have undergone a simi- 
lar mass accretion in the last 9-10 Gyr (z < 1 .5 - 2), to properly con- 
sider the observed variations of the central and the effective stellar 
mass densities of ETGs, one must take into account the following 
property: if the light profile of a galaxy and hence its stellar mass 
profile is represented by a Sersic profile (eq. 1) with index n > 2 
then the mass density within a fixed radius (e.g. 1 kpc) and the ef- 
fective stellar mass density are related exactly as shown by the data 
in Fig. 4, independently of redshift and of mass accretion experi- 
enced. This is shown in Fig. 4 by the red solid lines. The lower solid 



line represents the locus occupied on the [pe,pi] plane by galaxies 
described by a Sersic profile with index n = 4, having stellar mass 
10'" Mq and effective radii in the range 0.5 < R,, < 10 kpc. The 
middle and the upper solid lines represent galaxies with the same 
parameters but having a stellar mass 10" Mq and 4x10" Mq re- 
spectively. We see also that the lower the Sersic index the stronger 
the correlation between pi and pe as shown by the steeper dashed 
line obtained for an index n = 2 (normalized to 10" Mg). Both the 
stellar mass range and the effective radius range considered to com- 
pute the curves are those covered by our data. We see that the ef- 
fective stellar mass density of these simulated galaxies spans about 
three orders of magnitude while the central mass density pi is con- 
strained within an order of magnitude. Thus, any galaxy described 
by a set of parameters [n, A4,, within the above intervals (n > 2, 
M. = 0.1 - 4 X 10" Mo, R, = 0.5 - 10 kpc) will lie between the 
upper and the lower lines independently of redshift and of accre- 
tion history. In fact, our high-z ETGs fall between the two lines 
and spans the range bound by the two curves. In order to verify 
the independence of redshift of this evidence we used the ~900 
early-type galaxies at z < 0.05 selected from WINGS survey. We 
computed the central mass density pi and the effective mass den- 
sity pe of these ETGs using the same procedure described above 
and used for our high-z sample. The ~900 local ETGs are shown 
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in Fig. 4 as crosses. As expected, they fall in the region bounded 
by the two red lines covering the same distribution as the high-z 
sample. Fig. 4 also shows the 22 super-massive (0.3 - 5 x 10'^ Mq) 
ellipticals selected by Bemardi et al. (2006; 2008) on the basis of 
their high velocity dispersion (o-,, > 330 km/s) and studied by Tiret 
et al. (201 1) (open squares). It can be seen that these ETGs follow 
the same behavior shown by the other data: a central stellar mass 
density almost constant (but centered at a higher value according to 
the higher stellar mass range) and a large variation of p^. Finally, it 
is also important to note from Fig. 4 that the effective stellar mass 
density, the quantity which is expected to change dramatically (as 
the cube of the variation of R^) in the inside-out scenario from high- 
z to low-z, is not different in the two samples considered, in agree- 
ment with what is shown by our high-z sample in the lower panel of 
Fig. 5. This is more clearly shown in Fig. 6 where the distribution 
of the effective stellar mass density of our sample at 0.9 < z < 2 is 
compared to the distribution of p,, of the WINGS ETGs. The two 
distributions are equal as confirmed by the Kolmogorov-Smirnov 
test which provides a probability P=0.85. Thus, these results show 
two main facts. The first one is that, all the ETGs independently of 
redshift and of their accretion history have the mass density interior 
to a fixed radius within a narrow range of values if compared to the 
stellar mass density within the effective radius. This fact reflects a 
peculiar feature of the true intrinsic profile shape that they follow 
which is very well reproduced by a Sersic profile (see also Graham 
and Driver 2005). Since this behavior depends on the intrinsic pro- 
file shape of ETGs, it follows that also non-parametric estimators, 
e.g. Petrosian radius, would return qualitatively the same results. 
Thus, whether or not an ETG formed as a compact spheroid its 
central stellar mass density will fall within a narrow range of val- 
ues depending only on its total stellar mass. It does not depend on 
any past assembly history and has no connection with the possible 
inside-out growth. The other important evidence is that the effective 
stellar mass density of early type galaxies, the quantity expected to 
decrease dramatically in the inside-out growth scenario as the ef- 
fective radius increases, has not changed in the last 9-10 Gyr This 
implies that the large spread of and of R,, observed both in the 
local universe and up to z ~ 1.5-2, as shown in figures 2, 3 and 4, 
must originate earlier, at z > 2. 

4.2 Stellar mass densities and the evolution of ETGs at z > 2: 
toward tlie early-pliases of their assembly 

Figures 4 and 6 clearly show that the properties of the local ETGs 
in terms of stellar mass and stellar mass densities were already in 
place at z ~ 1.5-2. For this reason such properties must necessarily 
be the product of the evolutionary path followed by ETGs at z > 2, 
in the previous 2.5-3 Gyr i.e. during the assembly of most of their 
mass. Before making any hypothesis, let's see whether the general 
scheme of compact spheroids formation and evolution described 
above can justify the spread of the effective stellar mass density ob- 
served at z ~ 1.5-2. From figures 2 and 4 we see that, at a given 
mass, there are galaxies with effective radii (effective stellar mass 
densities) differing even by a factor 5-6 (more than hundred) even 
if all of them have nearly the same central mass density. Thus, in 
the hypothesis that ETGs formed first as compact spheroids, those 
galaxies that at 1 < z < 2 have large radii (5-6 times larger) with 
respect to others with equal mass, had to increase their size at an 
earlier epoch, at z > 2. Hence, the inside-out accretion should have 
been extremely efficient at z > 2, in the first 2-3 Gyr of the universe 
to produce so many large galaxies in so short time. Some authors 
have already discussed this point in different contexts (e.g. Nipoti 
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Figure 5. From top: the stellar mass M,, the effective radius Rp, the central 
stellar mass density pi and the effective stellar mass density pg of the 34 
ETGs of our sample at 0.9 < Zspec < 2 are plotted as a function of their 
redshift. No significant correlation is present. Pspearmen is the probability 
that the two quantities ai'e not con'elated as resulting from the Spearmen 
rank test. 
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Figure 6. The distribution of the effective steffar mass density of ETGs at 
0.9 < z < 2 (red histogram) is compared to the distribution of the effective 
stellar mass density of the local ETGs selected from the WINGS survey 
(blue histogram). The effective stellar mass density is the quantity expected 
to dramatically change in the inside-out growth scenario of ETGs. The two 
distributions are equal as confirmed by the Kolmogorov-Smimov test. 
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et al. 2009; Newman et al. 2011; Saracco et al. 2011) concluding 
that the required dry merging does not seem to be consistent with 
model predictions. Here, we can directly give an estimate of the 
merger rate needed on the basis of the observed data. Let us con- 
sider a galaxy of our sample having effective radius (effective stel- 
lar mass density) 4 (~60) times larger than another galaxy of equal 
mass at the same redshift. Let us assume that dry minor merger is 
very efficient in increasing the size of galaxies, as efficient as pro- 
posed by Naab et al. (2009) where the effective radius grows by the 
square of the change in mass instead of linearly. The larger galaxy 
has already increased its radius by a factor 4 during its life at z > 2. 
To increase the radius by a factor 4 the larger galaxy had to double 
at least its mass at z > 2. In the case of mergers involving galaxies 
with masses, e.g., M : m < 1 : 5 (1:10), the galaxy had to undergo 
five (ten) mergers in the previous ~ 2.5 Gyr (between 2 < z < 6) to 
double its mass. Thus, the resulting merger rate per galaxy would 
be about 2 (4) mergers per Gyr. If the final stellar mass at z ~ 1.5-2 
of the larger galaxy was 10" Mq the progenitor should have a mass 
M ~ 5x10'" Mq at z > 2. Using the merger rate calculator by Hop- 
kins et al. (2010) we estimated at z = 2.5 a merger rate per galaxy 
of about 6 X 10"' Gyr"' for masses in the range 5x10'" - 10" Mq. 
We considered dry merging, i.e. mergers involving gas fractions 
up to 10% of the total mass (gas -I- stars). The predicted merger rate 
would be even lower for lower masses (6 x 10""* Gyr"' for masses 
< 5 X 10'" Mo), for higher redshift (3 x lO""* Gyr"' at z = 3) and for 
lower mass ratios (4 x 10"' Gyr"' for M : < I : 10). Thus, the 
required merger rate at z > 2 as directly derived from the observed 
data to justify the data themselves in the hypothesis of accretion 
through minor mergers, is three orders of magnitude higher than 
the one predicted by models. Moreover, such high merger rates are 
excluded also from the merger rate derived from the observations 
of galaxy pairs at high-z. For instance, Lopez-Sanjuan et al. (201 1) 
estimate a minor merger rate R„,„, ~ 0.04 Gyr"' at z ~ 0.8 and they 
find that it should decrease with increasing redshift. Thus, at z ~ 1.5 
we should expect a minor merger rate R,„,„ < 0.04 Gyr"' , more than 
two orders of magnitude lower than the one required above. 

Some recent studies which find almost exclusively compact 
ETGs at high-z contrary to our observational evidence, consider the 
possibility that this difference is due to the higher redshift range 
covered by their samples. For instance, Cassata et al. (2010) and 
Szomoru et al. (2011) estimate the size of ETGs at 1.5 < z < 2.5 
using WFC3 data (0.13 arcsec/pix vs FWHM~ 0.15 arcsec) finding 
almost all compact ETGs. They select (and define) ETGs accord- 
ing to their colors and/or specific star formation rate. A possible 
suggested explanation for the higher fraction of normal ETGs ob- 
served at lower redshift (z ~ 1-1.5) is that some of the compact 
ETGs at high-z increase their size between z ~ 2.5 and z ~ 1.5, i.e 
in about 1.5 Gyr. However, if this was the case, following the same 
reasoning as above, a galaxy should have experienced from five to 
ten mergers in about 1.5 Gyr, requiring a merger rate of 3-7 Gyr"' 
at 1.5 < z < 2.5, even higher than the one excluded before. It seems 
then that the different fraction of compact galaxies they found with 
respect to our sample cannot be accounted for by this hypothesis 
but rather it should be attributed to other factors. For these works, 
particular attention should be paid to the different criteria used to 
define and hence to select early-type galaxies. In this regards, it is 
worthwhile to mention the work by van der Wei et al. (201 1). They 
investigated the stellar structure of massive, quiescent galaxies in 
the redshift range 1.5 < z < 2.5 using HST-WFC3 images. They 
found that their efi'ective radii were actually smaller than ETGs in 
the present-day universe with equal stellar mass. However, they es- 
timate that 65±15 per cent of the population of massive quiescent 



galaxies at z ~ 2 are disk dominated. Thus, the hypothesis of an 
inside-out growth driven by dry merging does not seem a viable 
mechanism to justify the large spread of the effective radius and 
of the stellar mass densities characterizing the population of ETGs 
at I < z < 2. To reinforce this conclusion it is also the detection 
of color gradients high-z ETGs (Gargiulo et al. 2011; Guo et al. 
2011) In particular, Gargiulo et al. (2012) detect color gradients in 
some but not all the ETGs of the sample at 1 < z < 2. Most impor- 
tantly, they find that the presence or the absence of color gradients 
is not related to the compactness of the galaxy: color gradients are 
presents both in compact and in normal ETGs as well as some of 
the ETGs, both compact and normal, do not show color gradients. 
Necessarily, inside-out growth driven by dry merging cannot be as- 
sumed as a general scheme of accretion but a different scenario 
must be hypothesized. We further analyzed the properties of high-z 
ETGs in order to gain new insights and constraints on the possible 
mechanism(s) able to produce the observed properties. 



5 CENTRAL STELLAR MASS DENSITY AND SCALING 
RELATIONS 

In the left panel of Fig. 7 the central stellar mass density pi is plot- 
ted versus the stellar mass M, of the galaxies. It is evident the tight 
correlation between the two quantities. The correlation is well re- 
produced by the relation 

pi = 2.35(±0.3) X 10' M"-^***""^'. (4) 

for our high-z ETGs (filled circles). We note that the local sample 
of ETGs (crosses) follow the same correlation with a best-fitting re- 
lation possibly slightly steeper (pi oc A4""). It is worth noting that 
Tiret et al. (201 1) found an opposite correlation using their data set 
(pi oc M'^"'^^), in contrast with our result and also with what the 
data of the other authors they used show (see their Fig. 6). There 
are two main differences between their work and the others which 
could justify this discrepancy. The first one is that they consider 
the sample of super-massive ellipticals (0.3 - 5 x 10'^ Mq) with 
velocity dispersion larger than 330 km s"' (Bernardi et al. 2006), 
while the other samples (as well as our own) include galaxies less 
massive with no cuts in velocity dispersion. The other difference is 
the stellar mass estimator used to derive the density within 1 kpc. 
Indeed, while we and the other authors have used the light profile 
to derive the fraction of stellar mass within a fixed radius, they use 
a dynamical mass computed by solving the Jeans equation. To di- 
rectly verify whether the different stellar mass range characterizing 
their sample can justify the different result, we derived the central 
stellar mass density of their ETGs according to our method (us- 
ing eq. 2) using the stellar mass they obtain from stellar population 
model (column 6 in their Tab. 2). Their sample is represented in 
figures 7 by open squares. It is evident in this case the agreement 
of their data with the others data in all the cases, showing that the 
different mass range and the selection in velocity dispersion are not 
the reasons of their opposite correlation. We believe that the use of 
the dynamical mass they define is the reason of the different result. 
Actually, their Fig. 4 shows the large scatter which exists between 
the dynamical mass and the stellar mass they derived and a relation 
between them slightly steeper than I : I which could be the reasons 
of the different behavior they find. The result obtained from this 
comparison shows that the correlation between the central stellar 
mass density and the stellar mass discussed above and the resulting 
considerations are valid for the population of ETGs independently 
of redshift, of mass and of velocity dispersion ranges considered. 
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Figure 7. Left - The central stellar mass density p\ is plotted versus the total stellar mass M, for our sample of high-z ETGs (large filled magenta and cyan 
circles) and for the WINGS sample of local ETGs (black crosses). These two quantities are tightly correlated independently of the redshift and the compactness 
of the galaxies. The correlation is well represented by eq. (4). Right - The fraction of the stellar mass within 1 kpc radius is plotted versus the total stellar 
mass of the galaxy. Symbols are the same as in the left panel. The existence of the correlation shown in the left panel and of the form of eq. (4) implies the 
correlation shown in this plot: the fraction of stellar inass contained within 1 kpc radius decreases with increasing total stellar mass of the system (eq. (5)). 



Because pi = 0.239A4i, it follows that also the central stellar 
mass increases with the total stellar mass of the galaxy according 
to the same scaling relation, Mi oc Ml'^^. Furthermore, from the 
above relation it follows that the fraction of the central stellar mass 
decreases with mass as 

Ml /M. =9.85^4;°"^ (5) 

Hence, if the stellar mass of the galaxy increases by a factor two 
then the stellar mass and the density within 1 kpc increase more 
slowly, by a factor 1.5. Consequently, the fraction of stellar mass 
contained within 1 kpc, MilMt, decreases by a factor 1.3 when the 
total mass double as shown in the right panel of Fig. 7. Finally, it is 
interesting to note that since the effective radius of ETGs increases 
with mass as R, oc M~°-^^ (e.g. Shen et al. 2003) and pi oc M;"-* 
(eq. 4) it follows that pi is nearly independent of R,,, as it is. This 
is another and independent way to show that pi and p^ are nearly 
uncorrelated, independent of redshift and of assembly history since 
the relations used to derive this correlation are valid both for high- 
z ETGs and for local ETGs independently of their compactness. 
Moreover, this indirectly shows that the small variation of the cen- 
tral stellar mass density of ETGs with respect to the effective mass 
density is an intrinsic characteristic of the Sersic's profile and has 
no regards with the inside-out growth. 

In all the cases considered above there are no differences be- 
tween the high-z sample and the local sample in spite of the ex- 
tra 9 Gyrs for local ETGs compared with high-z ETGs. All the 
ETGs follow the tight correlations shown above independently of 
their compactness and of their redshift. Hence, whatever evolution 
experienced by ETGs during the last 9-10 Gyr of their life, must 
be able to leave their properties unchanged. Such properties, al- 
ready in place at z ~ 1.5-2, must necessarily originate at earlier 
epoch, during the first 3 Gyr The large spread of the effective stel- 
lar mass density equally observed in the ETGs at z ~ 1.5 and in 
the local ETGs representing the large variation of R^ at fixed stellar 



mass, must originate in the early-phases of their formation. And the 
mechanism(s) responsible for mass growth of ETGs must proceed 
according to the above scaling relations. 



6 SUMMARY AND CONCLUSIONS 

We studied the relations between the stellar mass density within the 
effective radius p^ and the stellar mass density within a fixed radius 
of 1 kpc pi for a complete sample of 34 ETGs at 0.9 < Zspec < 2 
with stellar masses in the range 0.6 x 10'" - 3 x 10" Mq. As pre- 
viously found by other authors (e.g. Tiret et al. 201 1; Bezanson et 
al. 2009) we find that the central stellar mass density pi of ETGs 
varies just by a factor of two and it is similar to the one of local 
ETGs. This property of the high-z ETGs was supposed to be the 
sign of the spheroid formation scenario in which ETGs assemble 
at high-z as compact/dense spheroids and then add a low stellar 
mass density envelope growing efficiently in size and less in mass. 
The central stellar mass density and the stellar mass of the galaxy 
would remain nearly constant while the effective stellar mass den- 
sity would strongly decrease across the time as the effective radius 
increases. However, we find that the effective stellar mass density 
Pe of high-z ETGs varies by almost 3 orders of magnitude exactly 
as the local ETGs. We find indeed that high-z and low-z ETGs se- 
lected in the same stellar mass range populate the same locus in the 
[pe\Pi\ plane and that the distribution of the effective stellar mass 
density did not change since z ~ 1.5, i.e. in the last 9-10 Gyr This 
implies that the large spread of observed both in the local uni- 
verse and up to z ~ 1.5-2 must originate earlier, at z > 2. Moreover, 
we show that the small variation of the central stellar mass density 
with respect to the large variation of the effective density does not 
depend on redshift and/or on the past assembly history. On the con- 
trary, this feature is characteristic of the Sersic profile. Thus, the 
relation between the stellar mass densities of ETGs is not represen- 



© 2012 RAS, MNRAS OOO.fTlfTTI 



10 p. Saracco et al. 



Table 1. For each galaxy of the sample we report the spectroscopic redshift, the total apparent magnitude F850LPy;, resulting from the best fitting profile on 
the HST-ACS images, the eft'ective radius R,,, the Sersic index n, the effective stellar mass density p^, measured within the eft'ective radius, the central stellar 
mass density p\ measured within 1 kpc radius, the total stellar mass Al, and the stellar mass Mi within 1 kpc, the mean age of the stellar population. We 
remind that pi = 0.239A1i. The last column (Oil),,,,, reports those galaxies showing Oil emission. 
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tative and cannot be used to probe their past assembly history. It has 
no connection with the inside-out growth. We investigated whether 
the spread of the effective radius and of the effective stellar mass 
density observed in ETGs up to ~ 1.5 can be the sign of the inside- 
out growth taken place at higher redshift (z > 2). We find that the 
number of mergers needed to justify at z ~ 1.5 a spread of a factor 
4 in or about 60 in is from five to ten in the redshift range 
2 < z < 6, corresponding to a merger rate of about 2-4 Gyr"'. This 
merger rate is three orders of magnitude higher than the one pre- 
dicted by models and than the one derived from the observations of 
high-z galaxy pairs. 

As also discussed in Saracco et al. (2011), a scenario where the 
main starburst triggered by the dissipative gas-rich merging is the 
main driver of the stellar mass growth, could qualitatively account 
for the observations when supposing that gas cooling can modu- 
late the star formation. Indeed, assuming that dissipative gas-rich 
merger at high-z (z > 4 - 5) is the main mechanism of spheroids 
formation, compact ETGs would be a natural consequence of this 
mechanism when most of the gas at disposal is burned in the cen- 
tral starburst, that is when the gas is sufficiently cold to collapse 



toward the center and to ignite the main starburst. On the contrary, 
when the gas of the progenitors is not sufficiently cold and homoge- 
neous the resulting starburst would not be short, central and intense 
but longer and possibly composed of many episodes. The cooling 
time of the gas clouds and their orbital parameters could modulate 
the rate at which the starburst(s) are ignited, stoked and distributed 
across the galaxy. 

We find and define a tight correlation between the central stel- 
lar mass density of ETGs and the total stellar mass in the sense that 
the central stellar mass density increases with mass as Alr"^- We 
find that this correlation is general and it is valid for the whole pop- 
ulation of ETGs, independently of redshift, mass and/or of velocity 
dispersion ranges considered. Given this correlation it follows also 
that the fraction of the central stellar mass of ETGs decreases with 
mass. 

The large spread of the effective stellar mass density equally 
observed in the ETGs at z ~ 1.5 and in the local ETGs representing 
the large variation of at fixed stellar mass, must obviously origi- 
nates at z > 2. The fact that this property apparently did not change 
in the subsequent 9-10 Gyr as it is equally observed in the local 
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ETGs, suggest that it must have originated during the early phases 
of their formation. Moreover, whatever the evolution ETGs under- 
went during the last 9-10 Gyr of their life, such evolution must be 
able to leave unchanged their properties and the correlation found. 
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